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Exhibiting antiferromagnetic interlayer coupling in dilute magnetic semiconductor multilayers is essen- 
tial for the realisation of magnetoresistances analogous to giant magnetoresistance in metallic multilayer 
structures. In this work we use a mean-field theory of carrier induced ferromagnetism to explore possible 
(Ga,Mn)As based multilayer structures that might yield antiferromagnetic coupling. 
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1 Introduction 

The exciting new prospect of spin based electronics, known as spintronics, was initiated in 1988 with 
the discovery of giant magnetoresistance (GMR) in metallic multilayer structures 01 El- These structures 
consist of interposed ferromagnetic (FM) and non-FM layers. When adjacent FM layers are aligned in 
antiparallel directions, enhanced spin scattering of carriers causes an increased electrical resistance through 
the layers. Dilute magnetic semiconductors (DMS) such as (Ga,Mn)As are a novel class of FM materials 
which show many spintronic functionalities, and are considered promising candidates for future spintronic 
applications |3|. 

Antiferromagnetic (AFM) interlayer exchange coupling (IEC) in DMS based superlattices was theoreti- 
cally predicted in 1 999 using a k ■ p kinetic exchange model for carrier mediated ferromagnetism 1 4 1 . This 
approach considers delocalised charge and adds extra modulation induced by spin-polarised effects. A 
large magnetoresisticance (MR) was predicited due to the large difference in miniband dispersal between 
ferromagnetically and antiferromagnetically aligned layers. Recently, IEC has been further explored us- 
ing a tight binding model Q. This complementary microscopic approach, although not self-consistant, 
takes into account atomic orbitals for all the constituent atoms, leading to accurate descriptions of the band 
structure. Despite the different approaches used, both methods describe a Ruderman-Kittel-Kasuya-Yosida 
(RKKY) 2dkp interlayer coupling mechanism. 

In order to realise in DMS materials a phenomenon analagous to GMR in metals, but with potentially 
a much greater MR ratio, it is essential that AFM interlayer coupling is obtained. However, experimental 
work into (Ga,Mn)As based multilayer structures has only lead to reports of FM IEC |6|[7). This paper 
will describe part of a comprehensive review of the multidimensional parameter space available in DMS 
multilayer systems in order to identify optimal parameters for realising an antiferromagnetically coupled 
systems. 
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Fig. 1 The IEC E c for a multilayer where the magnetic layers are 2 monolayers thick (m = 2) and have a 2% Mn 
concentration. There is a uniform impurity concentration throughout the structure, (a) This is shown as a function of 
the average 3D carrier concentration, n^D, and the number of monolayers of non-magnetic layer, n. (b) Shown as a 
function of 2(dm+i)feF. The red curve is an estimate of the ideal RKKY range function. 



2 Theoretical Modelling 

The following calculations will be based on the mean-field theory of carrier induced ferromagnetism |4|. 
In this model the band structure is solved using a Kohn-Luttinger kinetic exchange Hamiltonian using a 
parabolic band k ■ p approximation with an additional term J pl j representing the p — d exchange interaction 
between Mn spins and hole spins. The limitation of this approach is that a single parabolic band approxi- 
mation is used, sacrificing full quantitative accuracy for qualitative descriptions of a wide range of systems. 
Subtleties of the band-structure and spin-orbit effects are neglected. However, qualitative agreement with 
the data published in 1 5 1 at least partially justifies this approach. 

In these calculations we use experimentally determined value of J p d = 55 meV nm 3 (5). A hole mass 
m* = 0.5m e and a spin of local Mn moments S = I at T = K. Thermodynamics is treated on a mean 
field level. This is done using the standard formalisation of the local-spin density approximation using the 
Kohn-Sham equations for inhomogeneous systems. 

For the calculated IEC energy, E c , positive values correspond to FM interlayer coupling to be energet- 
ically favourable, and negative values correspond to AFM interlayer coupling being the favoured configu- 
ration. 



3 Results 

3.1 GaAs spacer 

In the RKKY model of interlayer exchange the oscillations occurs as a function of dkF, where d is the 
separation between the two-dimensionsal magnetic planes and Uf is the Fermi wave vector jS). In our 
model we shall denote d n as the width of the non-magnetic layers, corresponding to d from the RKKY 
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model, and d m as the width of the magnetic layers. We shall also define the average Fermi wave vector kp 
as 



k F = (3^ 2 n 3D )3, (1) 

corresponding to the Fermi vector kp in the ideal RKKY model with a parabolic band. The average 
three-dimensional (3D) carrier concentration n 3 p> is defined as 



n 3D = - / n 3D (z)dz = - . (2) 

U-n+m J unit cell «n+m 

First we shall consider a multilayer structure close to the RKKY limit of infinitely thin magnetic layers 
surrounded by free unpolarised carriers. So, we shall use thin magnetic layers and a low magnetic moment 
concentration. In Fig. |l(a)| the IEC energy, E c , is plotted against the 3D carrier concentration, n^a, and 
the number of monolayers of GaAs in the non-magnetic spacer, n, where one monolayer has a thickness of 
0.283 nm. The magnetic (Ga,Mn)As layer is 2 monolayers thick and contains 2% Mn local moment doping. 
There is a uniform acceptor density throughout the structure which gives an average hole concentration of 
4.43 x 10 20 cm -3 . In this case there are oscillations as a function of both parameters, analogous to the 
dkp oscillations in the ideal quasi one-dimensional RKKY model. 

The data from Fig. |l(a)| is replotted in Fig. |l(b)| as a function of 2(d n +i)kp. Also plotted is the function 



where a is a scaling factor. This function is the asymptotic limit of the pseudo one-dimensional (ID) 
RKKY range function |8|. The strength of the interaction is expected to scale with the density of states, 
and in the ID case a ~ kp (10|. The different series of points on the graph correspond to the series of 
different n values from the original plot. For a given 2(d n +i)kp, the points with the largest magnitude 
are those with the greatest kp; this behaviour is consistent with the expected scaling of a with kp. The 
important point to note is the fact that, in order to have improved alignment of the curves, the oscillations 
were plotted using the parameter d n +i; the non-magnetic spacer thickness plus an additional monolayer. 
This qualitatively obtained correction factor is necessary because the magnetic layers are not longer 2D 
planes, so small spacers have proportionally less effect on the centre-to-centre distance of the magnetic 
layers. 



3.2 AlGaAs spacer 

To consider structures further from the the RKKY ideal, a band offset of 150 meV in the valence band is 
introduced. This represents a non-magnetic layer material of Alo.3Gao.7As and will cause a depletion of 
carriers from the non-magnetic layers, confining them to the magnetic layers. Fig. |2(a)| shows the IEC for a 
structure with a (Ga,Mn)As magnetic layer of 2 monolayers and an Alo.3Gao.7As non-magnetic layer with 
no acceptor doping. The peak FM and AFM coupling strengths are now stronger than with doped GaAs 
spacers seen in Fig. |l(a)| However, the 2d n kp oscillations are damped more rapidly than with the GaAs 
spacer, resulting in the second FM and AFM peaks being very weak. This additional damping occurs more 
rapidly with increasing carrier density, h 3 p>, than with with spacer thickness. As a result, the first AFM 
peak, which occurs at lower carrier concentrations for greater non-magnetic layer thinkness, is almost 
constant in value. This is in stark contrast to the GaAs barrier case, where the largest AF coupling when 
n = 10 is less than a quarter of the magnitude as when n = 2. This suggests there may be ways to retain 
strong AFM IEC at greater spacer thicknesses. However, in this case the carrier concentrations at which 
this occurs would be too low to expect (Ga,Mn)As to be FM. 
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With high magnetic layer thicknesses the RKKY type oscillations have almost completely disappeared. 
Fig. |2(b)| shows this for m — 8, with 2% Mn doping. In these cases oscillations occur almost exclusively 
with hole density, being almost independent of the spacer thickness. Furthermore these oscillations have 
lost their regular character and appear to have a beating type behaviour. This configuration offers some 
interesting possibilities. Since the oscillations are no longer strongly dependent on spacer thickness, the 
n,3D at which the peak AFM IEC occurs for a given d n is broadly constant. This is particularly interesting 
as a carrier concentration around 10 20 cm~ 3 is a technologically feasible quantity. It is worth noting 
that the magnetic layer thicknesses used in the experimental studies are many times thicker, at 25 and 50 
monolayers 1 6 1 and 20 nm (^70 monolayers) respectively 1 7 1 ; such large (Ga,Mn)As layers would not be 
expected to yield clear AFM IEC. 

4 Conclusion 

The composition and structure of (Ga,Mn)As based multilayers can have profound effects on the expected 
IEC. By examining possible compositions within the broad parameter space that these structures offer it 
is possible to identify different recipes for devices that might offer the possibility of demonstrating AFM 
interlayer coupling. 
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